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The theory of quantum feedback networks has recently been developed with the aim of showing 
how quantum input-output components may be connected together so as to control, stabilize or 
enhance the performance of one of the subcomponents. In this paper we show how the degree to 
which an idealized component (a degenerate parametric amplifier in the strong-coupling regime) 
can squeeze input fields may be enhanced by placing the component in-loop in a simple feedback 
mechanism involving a beam splitter. We study the spectral properties of output fields, placing 
particular emphasis on the elastic and inelastic components of the power density. 
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I. INTRODUCTION 

In the last two decades quantum physics has witnessed 
a remarkable convergence between theoretical models of 
interactions, particularly for open systems and measure- 
ment apparatuses, and experimental implementations of 
quantum engineering. The unifying framework has been 
to import and adapt the principles of control theory to 
the quantum domain. The advantages over traditional 
approaches show that quantum control will play a fun- 
damental role in emerging quantum technologies [TJ [2] . A 
variety of promising control techniques have been put for- 
ward [3] - (10] which extend open- loop paradigms (where 
control inputs are decided in advance) and measurement- 
based closed- loop paradigms (where feedback of observa- 
tions is used to determine the control inputs) . Real-time 
measurement-based feedback has been applied to adap- 
tive homodyne measurement |lll 112] to achieve measure- 
ment variances close to the standard quantum limit. 

Our interest lies in coherent quantum control, which 
is a non-measurement based feedback approach. Quan- 
tum feedback networks [13j [14] have emerged as a nat- 
ural class of objects with which to address assemblies of 
quantum input-output components so as to allow feed- 
forward and feedback connections. This offers a conve- 
nient framework to formulate problems in coherent quan- 
tum control and robust quantum control problems |15] - 
[21] . (We remark that the early formulation of coherent 
quantum feedback control due to Lloyd [22 deals with 
the direct interaction between system and its controller, 
as opposed to one mediated by quantum field processes. 
However, this may be treated as a special case of the 
network [T3].) 

An early application of feedback to enhance the squeez- 
ing of an (infrared) cavity mode was given by Wiseman et 
al. [23]. Here the mode is coupled to a second harmonic 
(green) mode which is subjected to a quantum nondemo- 
lition measurement. In contrast, we wish to examine the 
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FIG. 1: (Color online) Squeezing Device in a Feedback Loop 

squeezing of the input noise field by a cavity mode acting 
as an idealized squeezing device. Here the feedback is co- 
herent, rather than measurement-based, and we consider 
a set up involving a simple beam splitter to introduce 
the feedback loop. We shall work in the limit of instan- 
taneous feedback throughout. We shall be interested in 
the class of linear dynamical systems [23] , [T7] , [TH] , and 
indeed will study static components wherein the internal 
degrees of freedom have been eliminated. 

The degenerate parametric amplifier (DPA) is a well 
known non-linear device capable of squeezing input fields 
[25] -[27]. We follow the treatment of Gardiner [2SJ. For 
a single quantum input field coupled to a single cavity 
mode a with coupling strength ^fk and Hamiltonian 

tfDPA = |(a* 2 -a 2 ), (1) 

there is an approximate squeezing parameter given by, 
[25] section 7.2.9, 

r D PA = lnf— V (2) 
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Here the amplification is due to the specific choice of the 
Hamiltonian 7?dpa- 

Without feedback, the method of obtaining maximal 
squeezing for a degenerate parametric amplifier is to try 
and realize the Hamiltonian for the internal mode with 
parameter coefficient e as close to the threshold value 
(e = k) as possible, see [28 section 10.2. As originally 
noted by Yanagisawa and Kimura [23], the value of the 
effective damping for an in- loop mode, see FIG. [i] will 
depend on the reflectivity value a: 



k (a) 



1 



1 



(3) 



Our strategy is to use coherent feedback for a fixed degen- 
erate parametric amplifier (below threshold, and there- 
fore internally stable [18]) and tune the reflectivity of the 
beam splitter so as to select the degree of squeezing. 

The degenerate parametric amplifier is an idealized de- 
vice in which one assumes that k and e are large but with 
fixed ratio. We shall investigate the situation where both 
these parameters are finite. Also, we introduce additional 
quantum damping into the model to sec the effect of loss. 



II. QUANTUM FEEDBACK NETWORKS 
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FIG. 2: Input-output component 
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A. Components In Loop 

We may consider a feedback arrangement using a beam 
splitter, as in FIG. [3]below. The beam splitter is a static 
device which we take to be described by 



A single component consists of a quantum mechanical 
system, with Hilbert space f) driven by n quantum input 
processes b- m ,i (i = 1, • • • , ri), [55] [59], satisfying canonical 
commutation relations of the form [bi n i(t), b- ln jit')] = 0, 

[b in , i (t),b* Illj (t , )} = 8 ij s(t-t , ). (4) 

A schematic of a component appears in FIG. [2] 

The component is characterized by generator G — 
(S, L, H) where S = (Sij) is a unitary nxn matrix whose 
entries are operators on t) called the scattering coefficient 
matrix, L = (Li) is a column vector whose entries are op- 
erators on f) called the coupling coefficient vector, and H 
is a self-adjoint operator on f) giving the system Hamil- 
tonian. On the joint system-field space we have the uni- 
tary evolution process U (t) which satisfies the quantum 
Ito QSDE [30] 



dU (t) = I (Sy - Sij) dA in>ij (t) + LidB* n4 (t) 

- J2 LiSijdBiaj (t) -{\Y. + lH ^ dt \ U (*) ' 

id i ) 

(5) 

with U(0) = I. 

We encounter the integrated fields B in _i (t) = 
f Q b- m ,i{t')dt, B* nj (t) = B inj (t)* and A ia<ij (t) = 



( b ut \ _ T ( b in \ T _ ( a (3\ 

\ Win / V "out / ' " \(l 1/ J ' 

where T is taken to be a real-valued unitary matrix with 
determinant a — av — (3fj, = ±1. 

Suppose that G\ — (Sx,L\, Hf) is the generator of the 
(n = 1) component before the feedback connections are 
made. Once the component is in loop, in the limit of 
instantaneous feedback, we find an effective component 
with input b m and 6 ou t as indicated and generator given 
by Gi = (S\,Lx,Hi) with [13] 

Si = a + ^SlT 1 -^, (6) 
L x = (S^-vSlY 1 ^, (7) 
Hi = i/ 1 ° + Im{(L?)*(l-^)- 1 L?}. (8) 

The example above may be extended to include a loss 
mechanism describing coupling of the component to the 
environment, see FIG. [4] Prior to making connections 
we assume that the component is the four-port system 
(n = 2) with generator given by 

So = ( S Q l)' Lo= 

After feedback, the effective generator becomes G — 
(S,L,H) with 

S= ( S Q l) ' i= ' H = H i- 
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FIG. 3: Feedback using a beam-splitter 
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FIG. 4: Feedback with loss 



In the language of [12], [2] the effective generator is the 
concatenation G — (Si, L\, -ffi)ffl(l, L%, 0). We have rea- 
soned that since So is diagonal, there is no direct scat- 
tering between the inputs to the in loop device, and that 
concatenation of the effective lossless generator Gi with 
the loss mechanism G°. This however can be shown to be 
correct by utilizing the following construction from [13] : 
we note that the beam splitter itself can be understood as 
a static four-port component G' = (T, 0, 0) and the set up 
in FIG. [4] is then naturally identified as a Redheffer star- 
product arrangement of the two four-port devices, the 
effective generator for components in a Redheffer forma- 
tion is given in section 5.3 of |13j . and substitution into 
the expression gives precisely the generator G. 

We note that the relations we shall derive below for 
linear systems can be arrived at by algebraically elim- 
inating the internal fields vi n and w ou t- Whilst this is 
obviously easier than evoking the mathematical formula- 
tion of quantum feedback networks, we should point out 
that this is not entirely consistent and that the in-loop 
fields Vi n and v ou t are not canonical! Whilst this has 
been incorrectly interpreted elsewhere as a violation of 
the Heisenberg uncertainty relations, the reality is that 
the above description emerges from a regular model in 
which the commutation relations hold at all times how- 
ever the finite time delay of the feedback is taken into ac- 
count |13j . The in-loop fields are eliminated in the instan- 
taneous feedback limit and should then not be thought 



of as real physical fields. The algebraic arguments pre- 
sented here, however, reproduce the correct answer. 



III. LINEAR STATE-BASED INPUT-OUTPUT 
SYSTEMS 

We obtain a linear dynamical model in the case 
where our system is an assembly of quantum modes a a 
(a = 1, • • • , m) and the components of the generator take 
the special form Sij scalars, (using summation convention 
for repeated indices from now on) 



L 4 = CI 



C + n* 



H 



(9) 
(10) 



In this case it is possible to apply transform techniques to 
the dynamical equations. We define the transform fields 



/>oo 

6 [«] 4 / e - st b(t)dt. 
Jo 



Note that 



b* [s] = 



(/ e 
Jo 



s*t 



b(t)dty = b[s*y 



(ii) 



(12) 



Setting b- m ,i [s] — f Q °° e st h n i (t) dt, etc., we then obtain 
an input-output relation of the form 



6out,< [s]=Zr.( s )b in [ s } + E+( s )b* n Js} 



(13) 



where (s) are the transfer functions. (Here we ignore 
additional terms involving the system modes at initial 
time. This omission is justified when the model is stable.) 

It is convenient to introduce the doubled up nota- 
tion: for a vector x — (xi,--- ,xn) T , we write x = 
(xi, ■ • • Xpf, x\, • • • ,x* N ) T where T is transposition, for 

\ A B 

N x M matrices A, B we write A (A, B) = j ^ 

where U is entry- wise conjugation, [Ay] = [A*j . We also 
set A (A, B) — A (A^ , — _B T ) where f is the usual hermi- 
tian conjugation. We say that a matrix S = A (S-, 5+) 
is Bogoliubov, or symplectic, if it is invertible with 

S b = s- 1 . 

The transfer relation can be written as 

6 out [s] = E (s) b- m [s] 

with transfer matrix function 



E(s) 



E-(s) E+(s) 
E+(s)* S- (5)8 



Explicitly 
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S( S ) 



C [sh m - A) 



(14) 



where S = A (S, 0), C = A (C_, C+) where C T = [Cf a ] , 
A = -\C b C - iQ where Q = A (fl_,fi + ) with ft T = 



A. Analysis of the Spectrum 



In addition to the fields b [s] in ( 1 1 1 we also define past- 
field transforms 



*6 (t) dt 



(15) 



We remark that for input process as arguments, the fields 
b [s] and b* [s] commute with the fields c [s 1 ] and c* [a'] for 
all parameters s, s' since they involve integrals over future 
and past input fields respectively. 

The Fourier transform of a field b is then defined to be 



2tt ./_ 
1 



l b (t) dt 



b [0+ - iuj] 



The canonical commutation relations Q then imply 



that 



state we have 



(lu) , &i nj (a/) = SijS (lu — ui'). In the vacuum 



(6 in ,i [0+-iw] [0+-£o/]) 
(cin.i [0~ - «w] c* nj [CT - io;'] ) 



<%C+ (w + w'), 
(u; + lu') , 



where the Heitler functions are 



C+M 



or 



'dt, C- (") 



C± (lu) = n5(uj)±iPV 



° luJt dt, 
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In practice we shall only encounter the combination £ + + 
£_ = 2nd when calculating physical correlations, and not 
encounter the principle value contribution. In particular, 

(b in , i (uj)'bt n!j (uj')) = 6 ij S(u-u'), (16) 

have the sum of ^(&i n ,i [0 + — iou] b* n j [0 + + ioj']) 



and 



[o- 



;] c* n - [0 Likewise 



Ignoring the contribution from the initial value of the 
internal mode, the input-output relations for the past 



fields takes a similar form to (|T3|) namely 



-ij 0) Cin j- [s] + S+- (s) q* nJ [s] 



the only essential difference in the calculation being the 
sign change. Let us introduce the matrices 



Sin M = S,- • (-iw) , S± (lu) = a± (-iou) 



(17) 



then 



&out,i (w) = (u>) 6 inj - (cj) + 5+ (w) 6 in j- (-w)* . (18) 

We may therefore determine the correlation functions 
from the transfer functions given that the input is in the 
vacuum state: 



Kut.t H &out,j K)) = Mj (<*>) <5 (W - £«/) , 

feout,i(w)S outJ (w')) = Mij(u)(S(w + w')! (19) 



where 

Mj (w) = 5+ ( W )* 5+ (w) , M« (a;) = ( w ) 5+ (-a,) . 

(20) 

We note the straightforward identities 

Afij (cu)* = Afji (lu) , Mij (to)* = Mji (-cu) . 

We remark that 

3(-w) = A(5_(w),5 + H) 

and that this defines a Bogoliubov matrix for each real lu 
where it is well-defined, see [18] subsection V.C. In par- 
ticular, this ensures that the transformation from inputs 
to outputs is canonical, and the Fourier transform of the 
outputs satisfy a similar relation to |l6| ) . We see directly 
from (14 1 that lim^i^oo S (— iuS) = A(S, 0), or 



lim S- (lu) — S, lim 5+ (lu) 

— >oo |cj| — >oo 



0. 



(21) 



Definition: We say that a component is capable of 
spectral squeezing if the matrix M. (lu) is non zero for 
certain frequencies lu. In particular, given a vacuum in- 
put, we say that the ith mode is spectrally squeezed if 
■Ma (lu) 7^ for some lu. 

In the single input situation, the (lu) are complex- 
valued functions satisfying |«S_ (lu) | 2 — |<S+ (lu) \ 2 = 1. 
We then define the spectral squeezing function r (lu) by 
|<S_ (lu) I = coshr (lu), that is 

r H = l lD I^MI + I^MI 



ln{|S_ H| + |5 + H|}. 



(22) 
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B. Power Spectrum Density 



We define output quadratures by 

<Zou M (*> O) = e i9 b ut,i (t) + e- ie & out , 4 (t)* , (23) 

for fixed phases 9 £ [0,2tv). The integrated processes 
Qout,i(t,0) = f q ovi t,i(t',9)dt' are self-commuting for 
fixed 9 and different times t and indices i, and correspond 
to classical diffusion processes with Ito differentials sat- 
isfying 



dQ ut,i (t, 9) dQ ou tj (t, 9) — Sijdt. 
Following Barchielli and Gregoratti [32], we set 

Pij(u,e,T) = ~([ T e i ^q out , i (t 1 ,9)dt 1 
1 Jo 



(24) 



x / e 

o 

T 



— iwt2 



Qout,i (t2,9) dt 2 



e^^out.z (h,9) dh) 



</ e- wt2 q out , t (t 2 ,9)dt 2 ) 



Pl? el (u,9,T) = V 1] (lo,9,T)~V^lo 1 9,T), 
and define the power spectral density matrix to be 



lim P i:i (w,0,T), 



(25) 



whenever the limits exist, along with the clastic and 
inelastic components Pfj(uJ,9) = liniT^oo Pt}{u, 9, T), 



P™\lu, 9) = limr^oo P§°\u, 9, T) respectively. 



The Ito rule ( 24 1 implies that 



P%*\u,e) = Zih 



(26) 



and this may be interpreted by saying that the squeezing 
in the dynamic model comes entirely from the clastic 
component, and that there is no inelastic squeezing. 

The Fourier transform is then <j ut,i (w, 9) = 
e lS b ou t.i (uj) + e~ l9 b out ,i (— v)* and it is readily verified 
that, for vacuum input, 

<<Zou M (w, 0) Qout,j ((J, 6)) = Pa (w, 0) 8 (w + </) (27) 
where we obtain the explicit expression 

P lj {oj,9) = 5 ij +M ij (-u))+M ji {u) 

+ e 2l8 Mji (cj) + e~ 2ie Mij {-u})* ■ (28) 

C. Idealized Static Squeezing Components 

A static squeezing device is an idealized static compo- 
nent with input-output relation of the form (either in the 
time or transform domain) 



ij in J 



(29) 



where 



are constant coefficients such 



that S — A(S-, S+) is a Bogoliubov matrix. The out- 
puts are then a symplectic transformation of the inputs 
and therefore satisfy the canonical commutation rela- 
tions. 

In practice, such a device is realized approximately by a 
dynamical component, in a limiting regime. Specifically 
we would require in the Fourier domain that the coeffi- 



cients S^: (u) in (18 1 are approximately constant over a 



wide range of frequencies. 

Definition: We say that a sequence of models 
converges pointwise in transfer function if we have 
limfe-Kx, S fc (u) =E(u>). 

If the limit is a Bogoliubov matrix S = A (SL, S+) in- 
dependent of s, then we obtain a static device. In this 
case, if S+ = then «SL is unitary and the limit corre- 
sponds to a beam splitter with S = S-. The situation 
S+ can however arise as such limits, the DPA is an 
example, and we refer to such idealized components as 
static squeezing devices. This notion of convergence is 
weak since there is no quantum stochastic limit model 
for which we could obtain S+ =/= 0, specifically we would 
violate the requirement (211 common to all dynamical 



models considered up to this point. 

For the case of a single input {n = 1), the are 
scalars with the constraint |5_| 2 — \S+\ 2 = 1, which en- 
sures preservation of the canonical commutation rela- 
tions. The parameter r — cosh|5_| is referred to as 

the squeezing parameter. We then have | >S' j = coshr, 

15+1 = sinhr, and we find that the extremal squeezing 
ratios of quadratures by the device are e ±r . 

We should remark that the canonical transformation in 



equation ( 29 1 is a Bogoliubov transformation for a quan- 



tum field. There is a strict condition on when Bogoliubov 
transformations are unitarily implemented for infinite di- 
mensional systems (Shale's Theorem, [3T]) which are not 
met in this particular case. 



IV. THE DEGENERATE PARAMETRIC 
AMPLIFIER 

We now treat the specific example of a degenerate 
parametric amplifier. 



A. Lossy DPA, Open Loop 

We consider n = 2 input field processes driving a single 
(m = 1) mode with 



S 



1 
1 



L 



and H = i?DPA as in ([TJ. 
C+ = 0, fi_ = and fl+ 



Here C_ 



Therefore A 



6 



and we see that the system is Hur- 



k + 7 — e 
— e k + 7 

witz stable (that is, A has all eigenvalues in the negative 
half plane) if 



k + 7 > e. 



(30) 



We obtain the following expressions for Szp (s) 



P( S ) 



s 2 + 7s + 



2 2 2 
7 — K — E 

4 



2 2 2 
hi —7 — £ 

4 



~ + ( S ) 



where 



2P(s) 



K V K 7 
K7 7 



P( S )= s 



K + 7 



K + 7 



(31) 



(32) 



This gives the transfer function for the component in 
FIG. 3 prior to feedback connection. 



B. Lossy DPA, Closed Loop 



Let us take for definiteness the beam splitter matrix to 



be 



T{a) 



a (3 
(3 -a 



(33) 



where < a < 1 and (3 — y/1 — a 2 . Following our discus- 
sions in the previous section, the actual situation mod- 
eled in FIG. 3 is then given by replacing L\ by 



L\ (a) = P (1 — v) L\ = \J k (a)a 



(34) 



where n(a) = in accordance with (|3j. The trans- 

fer function is therefore of the same form as derived in 
( 3L 1 32 > but with k now replaced by k (a). 



C. Spectrum of the DPA Output 

We are interested in the input-output relation between 
&i nj i and 6 ut,i- Here, displaying the dependence on a, 
we compute 



r s 2 n(a) («(a) + 7) £n(a) 



AD(u),a) 



2L>(w,a) 



n(a) + 7 



(35) 



with 



D(w,a) = \P(-iu),a) \ 2 



K(a) + 7 + £ 



n{a) + 7 — £ 



(36) 



Note that in the lossless situation 7 = 0, we have 
(ui,a) 1 2 — \Sii (u>,a) \ 2 = 1 and therefore the identity 
|A4n(o7 a)| 2 = (J\fu(w,a) + l)jVn(w, a). In particular, 
we compute that the spectral squeezing function is 



D. The Static Limit of the DPA 

We consider a sequence of DPA models described by 
the parameters 7fc)fc>i with 



1 w 

r{io,a) = - In — 

Z 9 



k(q)+£ 

2 



(37) 



In this case we find that we find that the power spectral 
density V\\ (u, 9) is 



K fe = kit, £k = k£, 7^ = fe7, 

and consider the singular limit k — > 00. We note that 
(s) = S T (s//c) so the limit is equivalent to the low 
frequency limit. The limit transfer functions are inde- 
pendent of the transform variable s: 



£n(a) 



cos 29 S- = lim S W (s) 



D (uj, a) 

with the maximum squeezing at 9 - 
e 2r(a,ui) an( j mm j mum squeezing at 9 

p —2r(a,u) 



0, VufaO) 



(k + 7) 2 - £ 2 



K 2 -£ 2 



-2J~Rrj{n + 7) 



7 

-2J~KTf (k + 7) K 2 — 7 2 



5+ = lim S< fe) (s) = — 

fc->oo (t + 7) 



2e 



re V K 7 
fRrj 7 
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In particular, S — A (S-, S+) is a Bogoliubov matrix. 
The squeezing parameter for the limit Bogoliubov trans- 
formation is then 



In 



or equivalcntly r (0) in (22 1. This is of course the equa- 
tion 

The central issue here is that the asymptotic limit 
| a; | — > oo and the transfer function convergence limit 
k — > oo do not commute: denoting the spectral squeezing 
functions for the sequence of models as (lu) we have 
limfe^oo (uj) = r for all u, while lim^i^oo (uj) = 
for all k. 

We likewise find that the inelastic contribution to the 
power spectrum of the limit output quadratures is given 
by (9 = 0) 

7^(0 = 0)= r (K + 7 + £) \ x 
(k + 7 ) - £ 2 



(k + 7) 2 



+ e 2 - 2e( 7 



which has eigenvalues unity and 



(K + 7) 2 + e 2 + 2e( 7 -K) 
(38) 

' 1 and is there- 



[(«+7) 2 



fore positive definite as required. The matrix for 9 = ^ is 
obtained by replacing e by — e. We have seen that for the 
dynamic approximation we always have 'P™ 6 ' (oj, 9) = I, 
the identity matrix, for finite k. The limit situation on 
the contrary now has purely inelastic squeezing. 



V. FEEDBACK-ENHANCED SQUEEZING 

For an idealized static description of a lossless DPA 
when placed in loop as in FIG. 3, we find that the squeez- 
ing parameter is modified to 



= In 



k (a) + e 
k (a) — s 



with k (a) 



fr-^K. We observe that the critical value of 

l-\-a 

the reflectivity a is 



K 



ttcrit 



(39) 



Here k (a cr ; t ) = e, and the squeezing parameter diverges. 
The approximating dynamical model has spectral squeez- 
ing function 



r (lj, a„it) = x In ~ 

2 uj z 



(40) 



which possesses a logarithmic singularity at lu = for 
the critical situation. The open-loop system is stable if 
and only if k > e, while for the closed-loop system this 



is modified to k (a) > e. Therefore the infinite squeezing 
situation implies the onset of instability of the closed-loop 
amplifier. 

It is instructive to look at the lossy (7 > 0) closed- 
loop situation. The relevant description is then given by 
(351 and (|36[). Hurwitz stability requires n(a) + 7 > e. 



We assume that the dissipation is below the threshold 
value (7 < e) and that the open-loop system is stable 
(k + 7 > e), then the closed loop system is stable for a £ 
(0, a cr it), where the critical value is now 



«crit 



K — £ + 7 

K + e — 7 

Here the critical value solves k (a cr it) 
have 



Mn (w, a CT it) 



4 [lj 2 + e 2 ] lu 2 



e [e - 7] 



2 [1a 2 + e 2 ] , 



(41) 
7, and we 

(42) 
(43) 



which are finite for to ^ 0. Both expressions diverge in 
the limit k — > 00, however, when we replace the parame- 
ters by (Kfe,£fc,7fc). 



VI. CONCLUSION 

Coherent quantum feedback control offers an intrigu- 
ing potential to engineer physically interesting states, 
and achieve high performance for quantum devices. The 
feedback approach based on quantum measurement is 
limited to time scales set by the measurement apparatus, 
the computer estimating (filtering) the quantum state of 
the system based on the measurements, and the imple- 
mentation of the controls by the actuator based on the 
filtered state. In contrast, coherent control is limited by 
the time delays associated in light traversing the loop. 

We have shown that coherent feedback control can en- 
hance the capability of a device to squeeze quadratures by 
using an optical network involving a beam-splitter loop. 
Conversely, squeezing could be suppressed by altering the 
beam splitter, for instance, by reversing the sign of a, 
though there would arguably only enhancement would 
be desirable. 

By tuning the beam-splitter reflectivity, we can modify 
the effective damping of an in-loop degenerate parametric 
amplifier while leaving the amplifier The situation where 
the squeezing becomes infinite corresponds to the thresh- 
old value of the damping as discussed by Gardiner [28] , 
however, we observe that this 
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